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Abstract. Proposed silicon-based quantum-computer architectures have attracted attention because of their promise for 
scalability and their potential for synergetically utilizing the available resources associated with the existing Si technology 
infrastructure. Electronic and nuclear spins of shallow donors (e.g. phosphorus) in Si are ideally suited candidates for qubits 
in such proposals, where shallow donor exchange gates are frequently invoked to perform two-qubit operations. An important 
potential problem in this context is that intervalley interference originating from the degeneracy in the Si conduction-band 
edge causes fast oscillations in donor exchange coupling, which imposes significant constraints on the Si quantum-computer 
architecture. We discuss the theoretical origin of such oscillations. Considering two substitutional donors in Si, we present a 
systematic statistical study of the correlation between relative position distributions and the resulting exchange distributions. 



INTRODUCTION 

As semiconductor devices decrease in size, their physical 
properties tend to become increasingly sensitive to the 
actual configuration of dopant substitutional impurities 
[1]. A striking example is the proposal of donor-based 
silicon quantum computer (QC) by Kane [2], in which 
the monovalent ^'P impurities in Si are the fundamental 
quantum bits (qubits). This intriguing proposal has cre- 
ated considerable recent interest in revisiting all aspects 
of the donor impurity problem in silicon, particularly in 
the Si:^'P system. 

Two-qubit operations for the donor-based Si QC archi- 
tecture, which are required for a universal QC, involve 
precise control over electron-electron exchange[2, 3, 4] 
and electron-nucleus hyperfine interactions (for nuclear 
spin qubits). Such control can presumably be achieved 
by fabrication of donor arrays with accurate positioning 
and surface gates whose potential can be precisely con- 
trolled [5, 6, 7, 8]. However, electron exchange in bulk 
silicon has spatial oscillations [9] on the atomic scale due 
to the valley interference arising from the particular six- 
fold degeneracy of the bulk Si conduction band. These 
oscillations place heavy burdens on device fabrication 
and coherent control [10], because of the very high ac- 
curacy requirement for placing each donor inside the Si 
unit cell, and/or for controlling the external gate voltages. 

The potentially severe consequences of these prob- 
lems for exchange-based Si QC architecture motivated 
us and other researchers to perform further theoretical 



studies, going beyond some of the simplifying approxi- 
mations in the formalism adopted in Ref. [10], and incor- 
porating perturbation effects due to applied strain[l 1] or 
gate fields [12]. These studies, performed within the stan- 
dard Hei tier-London (HL) formalism [13], essentially re- 
confirm the originally reported difficulties regarding the 
sensitivity of the electron exchange coupling to donor 
positioning, indicating that these may not be completely 
overcome by applying strain or electric fields. The sen- 
sitivity of the calculated exchange coupling to donor rel- 
ative position originates from interference between the 
plane-wave parts of the six degenerate Bloch states as- 
sociated with the Si conduction-band minima. More re- 
cently [14] we have assessed the robustness of the HL 
approximation for the two-electron donor-pair states by 
relaxing the phase pinning at donor sites, which could 
in principle eliminate the oscillatory exchange behav- 
ior. Within this more general theoretical scheme, the 
floating-phase HL approach, our main conclusion is that, 
for all practical purposes, the previously adopted HL 
wavefunctions are robust, and the oscillatory behavior 
obtained in Refs. [10, 11, 12] cannot be taken as an arti- 
fact. 



SINGLE DONOR 

We describe the single donor electron ground state using 
effective mass theory. The bound donor electron Hamil- 




FIGURE 1. Electron probability density on the (001) plane 
of bulk Si for the ground state of a donor in Si within the 
envelope function approximation. The white dots give the in- 
plane atomic sites. The color scheme runs from red (high 
density) to purple (low density). The highest density is at the 
central site, corresponding to the impurity site. 



tonian for an impurity at site Rq is written as 

,JTo = jTsv + ^VO ■ 



(1) 



The first term, J^y, is the single-valley Kohn-Luttinger 
Hamiltonian [15], which includes the single particle ki- 
netic energy, the Si periodic potential, and the impurity 
screened Coulomb perturbation potential 



y(r) = - 



e|r-Ro 



(2) 



way the charge density at the donor site [where the 
donor perturbation potential Eq. (2) is more attractive] 
is maximum, thus minimizing the energy for i^r^ (r). 

In Eq. (3), (r — Rq) are envelope functions centered 
at Rq, for which we adopt the anisotropic Kohn-Luttinger 
form, e.g., for jJ. = z, F^i^) = exp{— [(x^ +y^)/a^ + 

/b^]^^'^}/ Vko^. The effective Bohr radii a and b 
are variational parameters chosen to minimize Egy = 
{WRo\'-^^sv\WRo)^ leading to a = 25 A, = 14 A and 
Esv ^ —30 meV when recently measured effective mass 
values are used in the minimization [10]. The periodic 
part of each Bloch function is pinned to the lattice, 
independent of the donor site. 

The J^/^sv ground state is six-fold degenerate. This de- 
generacy is lifted by the valley-orbit interactions [16], 
which are included here in J^o, leading to the nonde- 
generate (A i -symmetry) ground state in (3). Fig. 1 gives 
the charge density [^/({^(r)^ for this state, where the pe- 
riodic part of the conduction band edge Bloch functions 
were obtained from ab-initio calculations, as described 
in Ref. [14]. The impurity site Rq, corresponding to the 
higher charge density, is at the center of the frame. It is 
interesting that, except for this central site, regions of 
high charge concentration and atomic sites do not nec- 
essarily coincide, because the charge distribution period- 
icity imposed by the plane-wave part of the Bloch func- 
tions is iTt/kjx, incommensurate with the lattice period. 

The oscillatory behavior of the single donor wave 
functions in Si, illustrated in Fig. 1, is well established 
experimentally [17] and theoretically [14, 18]. This be- 
havior does not bring significant consequences for con- 
ventional applications in Si-based devices (n-doped Si). 
A recent study of the single-qubit operations (A-gate) in 
the Kane QC shows that the A-gate operations do not 
present additional complications due to the Si band struc- 
ture interference effects [19]. 



For shallow donors in Si, we use the static dielectric con- 
stant £ = 12.1. The second term of Eq. (1), J^o^ rep- 
resents the inter-valley scattering effects due to the pres- 
ence of the impurity which breaks the bulk translational 
symmetry. 

The electron eigenfunctions are written on the basis 
of the six unperturbed Si band edge Bloch states 0^ = 
M^(r)e''''' '' [recall that the conduction band of bulk Si 
has six degenerate minima {jj. = 1, ... ,6), located along 
the F— X axes of the Brillouin zone at |k^| ~ 0.85(2;r/a) 
from the F point] : 



var„ (r) = -= x; (r - RoK (i-y 
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(3) 



The phases of the plane-wave part of all band edge Bloch 
states are naturally chosen to be pinned at Rq: In this 



DONOR PAIR EXCHANGE COUPLING 

The HL approximation is a reliable scheme for the well- 
separated donor pair problem (interdonor distance much 
larger than the donor Bohr radii) [13]. Within HL, the 
lowest energy singlet and triplet wavefunctions for two 
electrons bound to a donor pair at sites R^ and Rg, are 
written as properly symmetrized combinations of i/ra 
and \//Rg [as defined in Eq.(3)] 
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where S is the overlap integral and the upper (lower) 
sign corresponds to the singlet (triplet) state. The energy 
expectation values for these states, £/ = 
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FIGURE 2. Exchange coupling between two phosphorus 
donors in Si along the indicated directions in the diamond struc- 
ture. Values appropriate for impurities at substitutional sites are 
given by the squares. The dashed line in the R\\ [100] frame is 
a guide to the eye, indicating that the oscillatory behavior may 
be ignored for donors positioned exactly along this axis. 



gives the exchange splitting through their difference, 
J — Et — Es- We have previously derived the expression 
for the donor electron exchange splitting [11, 14], which 
we reproduce here: 

AR) = 4LAv(R)cos(k^-kv)-R, (5) 

where R = Ra — Rb is the interdonor position vector 

and ^|jv(R) are kernels determined by the envelopes 
and are slowly varying [10, 11]. Note that equation (5) 
does not involve any oscillatory contribution from (r), 
the periodic part of the Bloch functions [12, 14]. The 
physical reason for that is clear from (3): While the 
plane-wave phases of the Bloch functions are pinned to 
the donor sites, leading to the cossine factors in (5), the 
periodic functions are pinned to the lattice, regardless 
of the donor location. 

The exchange energy calculated from Eq. (5) for a pair 
of donors as a function of their relative position along the 
[100] and [110] crystal axis is given in Fig. 2. This figure 
vividly illustrates both the anisotropic and the oscillatory 
behavior of J{R), which is well established from previ- 
ous studies [9, 10, 11, 12]. It is interesting to note that 
for substitutional donors with interdonor position vectors 
exactly aligned with the [100] crystal axis, the oscilla- 
tory behavior may be ignored in practice, as indicated 
by the dashed line in the figure. This behavior is quali- 
tatively similar to the exchange versus donor separation 



dependence assumed in Kane's proposal [2], where the 
Herring and Flicker expression [20], originally derived 
for H atoms, was adapted for donors in Si. Therefore 
one might expect that reUable exchange gates operation 
would be possible if donors are exactly aUgned along the 
[100] crystal axis. 



NANOFABRICATION ASPECTS 

Aiming at the fabrication of a P donor array accurately 
positioned along the [100] axis, and given the current de- 
gree of control in substitutional P positioning in Si of a 
few nm [5, 6, 7, 8], we investigate the consequences of 
interdonor positioning uncertainties in the values of the 
corresponding pairwise exchange coupling. We define 
the target interdonor position Rj along [100], with an ar- 
bitrarily chosen length of 20 lattice parameters (~ 108.6 
A) indicated by the arrow in Fig. 2. The distributions for 
the interdonor distances /? = |Ryi — Rb| when R^ is fixed 
and Rb "visits" all of the diamond lattice sites within a 
sphere centered at the target position are given in Fig. 3. 
Different frames give results for different uncertainty 
radii, and, as expected, increasing the uncertainty radius 
results in a broader distribution around the target dis- 
tance. Note that the geometry of the lattice implies that 
the distribution is always centered and peaked around Rt, 
as indicated by the arrows. The additional peaks in the 
distribution reveal the discrete nature of the Si lattice. 

The respective distributions of exchange coupling be- 
tween the same donor pairs in each ensemble is presented 
in Fig. 4, where the arrows give the exchange value at the 
target relative position: 7(Rr) ~ 0.29 meV. The results 
here are quahtatively different from the distance distribu- 
tions in Fig. 3, since they are neither centered nor peaked 
at the target exchange value. Even for the smallest uncer- 
tainty radius of 1 nm in (a), the exchange distribution is 
peaked around 7-^0, bearing no semblance to the inter- 
donor distance distributions. Increasing the uncertainty 
radius leads to a wider range of exchange values, with a 
more pronounced peak around the lowest J values. 

From the perspective of current QC fabrication efforts, 
~ 1 nm accuracy in single P atom positioning has been 
recently demonstrated [6], representing a major step to- 
wards the goal of obtaining a regular donor array embed- 
ded in single crystal Si. Distances and exchange coupling 
distributions consistent with such accuracy are presented 
in Figs. 3(a) and 4(a) respectively. The present calcula- 
tions indicate that such deviations in the relative posi- 
tion of donor pairs with respect to perfectly aligned sub- 
stitutional sites along [100] lead to order-of-magnitude 
changes in the exchange coupUng, favoring 7 ~ val- 
ues. Severe limitations in controlling J would come from 
"hops" into different substitutional lattice sites. There- 
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FIGURE 3. Interdonor distance distributions for a target rel- 
ative position of 20 lattice parameters along [100] (see arrows). 
The first donor is fixed and the second one "visits" all of the Si 
substitutional lattice sites within a sphere centered at the target 
position, with uncertainty radii (a) 10 A, (b)20 A and (c)50 A. 
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FIGURE 4. Exchange distributions for the same relative po- 
sition ensembles in Fig. 3. The arrow indicates the target situ- 
ation. Contrary to the distance distributions, the exchange dis- 
tributions are not centered or peaked around the target value. 



fore, precisely controlling exchange gates in Si remains 
an open challenge. 
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